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ABSTRACT 

We develop a new fast-diffusion approximation for the kinetics of deposition of 
extended objects on a linear substrate, accompanied by diffusional relaxation. This 
new approximation plays the role of the mean-field theory for such processes and is 
valid over a significantly larger range than an earlier variant, which was based on a 
mapping to chemical reactions. In particular, continuum-limit off-lattice deposition is 
described naturally within our approximation. The criteria for the applicability of the 
mean-field theory are derived. While deposition of dimers, and marginally, trimers, is 
affected by fluctuations, we find that the A;-mer deposition kinetics is asymptotically 
mean-field like for all A; = 4, 5, ... , oo, where the limit k — )• oo, when properly defined, 
describes deposition-diffusion kinetics in the continuum. 
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1. INTRODUCTION AND DEFINITIONS 

Recent experimental advances [1,2] in the studies of kinetics of formation of amor- 
phous layers of monodispersed submicron colloid particles and proteins on fiat surfaces 
have stimulated much theoretical interest in models of irreversible growth [2] . Mono- 
layer growth by random sequential adsorption was the earliest such system to be 
considered [2,3], and many theoretical results have been accumulated over the years. 
More recent efforts have been focused on multilayer structures, on results which can 
be obtained by modern computer simulations, and on allowing for relaxation in the 
deposit [4,5]. The subject of this work is diffusional relaxation in monolayer deposition. 

In colloid deposition on surfaces, relaxation is typically slow. Diffusional relax- 
ation can, however, be observed experimentally in protein monolayers and multilayers 
[6]. Furthermore, multiunit-long, "sliding" molecules can be attached to DNA suggest- 
ing consideration of one- dimensional lattice and continuum models. Related models of 
deposition of the epitaxial type (i.e., crystalline-structure forming), accompanied by 
diffusional and other relaxation processes, were considered in the literature [7]. 

Theoretical study of the irreversible adsorption with diffusion encounters several 
difficulties. Firstly, numerical simulations have proved exceedingly resource- consuming 
already even for one- dimensional models [4,5]. Secondly, even in the absence of de- 
position, systems of diffusing hard-core extended objects are not fully understood in 
dimensions higher than one, and are the subject of much current research and contro- 
versy [8]. Indeed, both the dynamics and equilibrium state of hard-core systems show 
complicated cooperative effects such as slow "glassy" relaxation, phase transitions and 
phase separation, etc. 

Results of numerical simulations of lattice models in Id were interpreted [5] within 
a certain mean-field type framework which will be described shortly. First, we have 
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to introduce some notation. Consider deposition of extended objects of size £ on a 
line, with the attempt rate R (per unit length and time). The attempts are uniformly 
distributed, but only those are successful which fall in the space free from objects 
deposited earlier. Initially, the substrate [Id line) is empty. The deposited objects hop 
randomly, with the single-object (dilute-limit) diffusion constant D. Lattice models of 
A;-mer deposition are defined by further introducing the space mesh in steps of 

b = l/k, (1.1) 

and requiring that both the deposition events and the hopping keep the extended 
objects aligned and registered with the resulting Id lattice. For instance, the deposition 
attempt rate is Rb per lattice site (per unit time). 

The limit D = corresponds to deposition without diffusion, and provided k > 
1, the coverage (fraction of length covered), 0[t), approaches a certain A;-dependent 
jamming value ^(oo) < 1 for large times, t. This Id deposition model is in fact exactly 
solvable [3]. The jammed state is formed because there are gaps of lengths 6, . . ., 
(A; — 1)6, which cannot be filled by deposition. There exists a well defined continuum 
limit A; — )• oo, for all quantities of interest. 

However, for nonzero diffusion rate the coverage eventually reaches 100%, for large 
times, because diffusion allows the formation of large gaps within which a fraction of 
deposition attempts can succeed. Exact solution is no longer possible, even in Id. 
Instead, the following line of argument has been advanced [4,5]. For large times, with 
6 near 1, most of the free length will be in isolated single-site gaps of length b. These 
can be brought together by diffusion to form large empty gaps, of at least k lattice 
spacings, which can with some finite probability be covered by a depositing object (or 
be broken up again by diffusion). The asymptotics of the difference (1 — ^) for large 
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times can then be related to the decrease in density of the diffusion-hmited chemical 
reaction 



kA inert , (1.2) 

where A are diffusing reactants on a line, forming inert species, i.e., effectively annihi- 
lating, according to (1.2) with some finite probability on each encounter of k reactants. 

Such reactions have been studied in the literature [9]. The appropriate results 
will be reviewed later. For now it suffices to point out that for A; = 4,5, . . ., and with 
some reservations (to be addressed later) for A; = 3, one can use the mean-field rate 
equation which in terms of the coverage takes the form 



- c:, j,R£{l - 6)' . (1.3) 

While this approximate power-law behavior was confirmed by numerical simulations 
[5] for k = 3,4,5,10, there is an important limitation to the rate equation (1.3): it has 
no obvious continuum limit as A; — )• oo. 

Indeed, no matter what A;-dependence one conjectures for the dimensionless effec- 
tive rate constant 7^, the right-hand side of (1.3) cannot approach a definite function 
of ^ as A; — )• 00. This difficulty was also suggested by a preliminary observation based 
on the data for A; = 3, 4, 5, 10, that the times t from which (1.3) applies increase rapidly 
with k. 

The main objective of the present work is to repair this deficiency of the mean- 
field theory in Id. Thus, in Section 2 we introduce a new approximation scheme, still 
mean-field in nature, which allows the derivation of improved rate equations which are 
valid in a significantly larger regime and which apply in all the limiting cases. Various 
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implications of this new mean-field description will be detailed in Sections 3 and 4, 
with the latter devoted to the continuum limit. The original rate equation (1.3) turns 
out to hold only for 

t:>[{k-i)m]-^e''-\ (1.4) 

which explains its failure to describe the continuum limit A; — )• oo. 

As with any mean-field approximation, fluctuations are expected to dominate the 
large-time behavior in some cases, such as for k = 2. The appropriate criteria are 
derived and discussed in Section 5. Finally, Section 6 is devoted to summary and 
discussion of those new aspects that must be considered in order to extend the theory 
to d > 1, which, however, is not attempted here; see earlier remarks in this section 
and the discussion in Section 6. 

2. FAST DIFFUSION AND DIFFUSER SPACE 

It is convenient to relate our system to a lattice model of hard-core particles 
(monomers) of length single lattice spacing 6, each attempting to hop with rate H = 
AD/b^ per unit time to the left or right neighbor lattice site (rates H/2). We will 
use the term "diffusers" to refer to these particles. Hopping attempts succeed only if 
the hard-core constraint is not violated. The precise dynamics of diffusers is in fact 
not of interest to us here. The only property we will use is that for large times the 
system of hard-core particles approaches a steady state in which they are distributed 
statistically uniformly along the line [10]. 

Consider now a large length L in our original deposition model. It is convenient to 
take L to be the system size, in this section. We disregard finite-size effects by implying 
the thermodynamic limit first, in all the "intensive" expressions. The length L contains 
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OLJI extended objects, each k lattice spacings long. Disregarding deposition, diffusion 
in this system is the same as in a single-site-diffuser system in which all the extended 
objects of length I = kh are replaced by diffuser particles of length h. The total length 
available to these diffusers is 

K = {i-e)L + eLik= (^-e + ^e^L . (2.1) 

The density of diffusers (per site of the diffuser-space lattice) is thus 



' = TJh=i-e + \e ■ ^'-'^ 

The gap size between every pair of neighboring diffusers can be mfe, where m = 
0,1,.... Note that m = corresponds to contact, while the upper limit is taken 
to be infinity due to our disregard of finite-Z effects. Diffusion makes particle loca- 
tions uncorrelated. In terms of the gap sizes, the appropriate statement is that the 
(normalized) probability to find gap size m is 



Prob(m) = p(l -p)™ . (2.3) 

Of course, the result (2.3) applies only for large times. However, we are proposing 
to use (2.3) as the approximate, fast-diffusion, mean-field form of the gap size distribu- 
tion for the actual deposition model. The implications of this approximate assumption 
will be worked out in this and the next two sections. We will then proceed to derive 
the conditions for its asymptotic validity, in Section 5. Let us point out, however, that 
in the opposite extreme of no diffusion at all, exact results for the gap size distribution 
can be derived [3] although they were analyzed in detail only in the continuum limit 
A; — )• oo. These results will be described in Section 4. 
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We now return to the original full-line problem, and assume that (approximately) 
the gap sizes are distributed according to (2.3), where the time-dependence enters via 
the ^-dependence implied by (2.2). The total number of gaps (including m = 0) is 
equal to the number of objects. However, only gaps of size m > k offer (m — A; + 1) 
distinct locations for deposition. Since the deposition attempt rate per lattice site is 
Rb, the overall average rate of successful depositions in Z, per unit time, is 

u; = {Rb){eL/i) f^{m-k + l)p{l - p)™ = ^^^^^ ~ . (2.4) 

m=k 

In each deposition, the coverage is incremented by i/L. Thus, we obtain the relation 

d0 IR0{1 - pf 



dt kp 

Finally, we use (2.2) and introduce the dimensionless time variable 



(2.5) 



T = R£t, (2.6) 
to write the equation for the coverage in the fast-diffusion approximation, 

<"> (1 - (27) 

3. LATTICE DEPOSITION WITH FAST DIFFUSION 

In this section we consider implications of the relation (2.7) for fixed A;, i.e., for 
lattice deposition. Some conclusions are valid also in the limit A; — )• oo which, however, 
will be considered separately in the next section. We assume k = 2,3,4, .. . since the 
case A; = 1 is trivial. In fact, (2.7) is exact for A; = 1. 
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The curve 0[T) obtained by integrating (2.7) with the initial condition ^(0) = 
provides the upper bound (corresponding to D = oo) for any deposition process with 
fixed diffusion rate (recall that the deposition rate was absorbed in the time scale); 
see (2.6). The exact solution for Z) = 0, available as quadrature [3], provides the lower 
bound. The two curves are in fact quite close for small and intermediate times as 
illustrated in Figure 1 for A; = 5. 

For large times, however, the behavior is quite different. The D = curve ap- 
proaches the jamming value < 1, at the exponential rate, A;~-'e~-^/*', see [11]. However, 
the fast-diffusion result approaches the full 100% coverage, according to 

1 - e{T) ~ A;-^ [(A; - . (3.1) 

Interestingly, this power-law behavior sets in, and (2.7) can be replaced by the "reac- 
tion" type rate equation (1.3) with jk = k^~^ ■, for times 

T > (A; - l)-^e*'-\ (3.2) 

which become quite large as k increases. This time scale is essentially determined by 
the condition (1 — 0) <^ imposed to have no variation due to the denominator in 
(2.7), for 6 near 1. However, from (2.2) one can easily check that the same condition is 
also obtained from [1 — p) <^ 1, within the fast-diffusion formulation. The latter condi- 
tion ensures single-lattice-spacing gap dominance, as discussed in Section 1, provided 
that the gap size distribution is approximated by (2.3). 

The full result (2.7) is, however, more general than the rate equation (1.3); its 
precise limits of applicability will be derived in Section 5. As found for chemical 
reactions [9], and more generally for critical phenomena, the mean-field results can be 
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Fig. 1. The A; = 5 curves 0[T) for D = oo (upper) and D = (lower). The broken 
hne marks the hmiting, jamming value of the D = coverage for large times. The 
D = oo coverage curve approaches 1 as T — )• oo; see Section 3. 

used in those cases when local fluctuations do not dominate the dynamics. 

However, even for the description of dynamics dominated by the uniform fluc- 
tuations describable asymptotically by effective-field type approaches some caution is 
called for. Indeed, one can claim that relation (2.7) provides the correct functional 
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form of the asymptotic large time behavior, with the "large times" defined much less 
restrictively than when using (1.3); see Section 5. However, for any finite D (and 
A; > 2) to actually fit the data to the mean-field form, we must use it with "effective," 
renormalized rather than "bare" parameters. In (2.7) the only parameter was the 
overall rate. Phenomenological fits of data obtained for Z) < oo, should be done with 
the right-hand side multiplied by an additional adjustable rate parameter of order 1, 
before integrating (2.7). This parameter will approach 1, and the time range for which 
(2.7) applies will extend down to T = 0, only in the limit D — )• oo. 

Due to the limited numerical data available in the literature on the deposition 
problem [5] we have not attempted systematic tests of the above expectations. Some 
preliminary checks were made for k = 10. However, in a forthcoming publication [12] 
we report a similar improved mean-field approach formulated for the chemical reactions 
(1.2), as well as detailed numerical studies of the applicability of such "fast-diffusion" 
mean-field theories in Id. 



4. CONTINUUM DEPOSITION 

The A; — )• oo limit of (2.7) is easy to work out. 



d0_ 
dT 



(l-^)exp(^-^) . (4.1) 



Thus the fast-diffusion approximation is well-defined. The approach to the limiting 
coverage is quite slow: 



This behavior is quite different from the other extreme of no diffusion at all. The 



coverage in the latter case reaches the jamming value (< 1) at the rate 




see [11] 



for details. 



The diffuser space concept was useful in our derivation, Section 2. However, it 
cannot be used directly in the continuum because the lattice spacing b vanishes. It 
is useful to describe how the gap size probability in the fast diffusion approximation 
can be analyzed directly in the original space. We will also compare the fast-diffusion 
[D = oo) expression to the result in the D = case. (The main conclusions also apply 
for A; < oo. As already mentioned, the discussion is presented in this section because 
the D = gap size distribution has been studied in detail only in the continuum limit 
[3,13].) The diffuser space is nevertheless an important tool in the derivation of the 
limits of validity of the mean-field approximation; see Section 5. 

To analyze the gap sizes in continuum directly, in the fast-diffusion approxima- 
tion, we note that for uncorrelated object locations, the gaps must be exponentially 
distributed. (As we just pointed out a similar line of argument can be carried out 
in the discrete case.) Let g > denote the gap size (formerly mfe), and G{g) the 
density of gaps of size g to g -\- dg per unit length of the system, so that G has units 
of l/(length)^. The decay width and the coefficient of the exponential can be fixed 
from the condition that the total number of gaps be equal that of objects, which can 
be reduced to 




(4.3) 



and another condition: that the total length of all gaps be fraction (1 — ^) of the system 



length 




(4.4) 
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The resulting exponential distribution function, 

can be also shown to represent the correct A; — )• oo limit of the gap probability derived 
in Section 2 by "diffuser" arguments. 

For r — )• oo the exponential (4.5) collapses to the point g = This is in contrast 
to the Z) = distribution which only develops [3,13] a weak (logarithmic) integrable 
singularity at = 0, and is finite up to = I. For g > I ii is zero, which is just 
the manifestation of jamming: large gaps are filled up but the small ones are left 
over. For finite times, the singularity at the origin is not sharp [3,13]. There is a 
discontinuity in slope aX g = I, and in fact for g > I the functional form is exponential 
in g. Specifically, for large times, on approach to jamming, the Z) = distribution 
follows ~ exp [—{T I l)g], whereas the fast-diffusion exponential (4.5) behaves for large 
times as ~ exp [—{XnT I l)g]. 

The main effect of diffusion is therefore to increase the weight of large gaps at 
the expense of small gaps, which is particularly effective in the late stages of the 
deposition. These considerations actually suggest possible directions of improving the 
fast-diffusion approximation. A natural theoretical development, not attempted here, 
would be to write kinetic equations for the gap size distribution, incorporating the 
time- variation due to both diffusion and deposition. 

5. CRITERIA FOR APPLICABILITY OF THE MEAN-FIELD THEORY 

In order for the mean-field type approximations to hold, the system behavior 
must be dominated by fluctuations in various overall averaged quantities. If instead 
local density fluctuations dominate the deposition process, then mean-field theory will 



(1 - 6)1 



(4.5) 



not be applicable. Thus, for mean-field theory to apply, density fluctuations due to 
successful deposition attempts must be effectively smoothed out by diffusion, over a 
certain characteristic length scale L (which in this section no longer denotes the system 
size). 

For short times the substrate is nearly empty, the deposited objects are uncorre- 
lated and so L will be small. We expect L to become large for large times. In order 
to estimate Z, we make two observations. Firstly, the average distance between two 
neighboring empty sites is 



This distance becomes large for 6 near 1. Secondly, each successful deposition elim- 
inates k empty sites, which were at that instant part of a gap of m > A; sites. The 
overall gap distribution will thus be appreciably "perturbed" over a length 



L = 0{k)AL , (5.2) 

which represents the average size of the region which would have been occupied by 
0{k) empty sites after diffusive spread, had there been no deposition event. 

In order to have a more definite expression, we note that for A; = 1 there are no 
correlations at all, and formally, L = 0. Thus we replace the 0{k) factor by (A; — 1), 

This prefactor choice is, however, purely phenomenological. 

For any region of length Z, the rate u; of successful deposition events in that 



region is given by (2.4). The time scale it takes the diffusion to equihbrate density 
perturbations in length L is 



(5.4) 



where A is the diffuser-space equivalent of Z, see (2.1). 



The condition to have effective diffusional equilibration is 



a;T < 1 



(5.5) 



which upon collecting all the relations and definitions introduced above reduces after 
some algebra to 



The right-hand side of (5.6) is shown in Figure 2 for k = 2,3,4,6 and oo, where 
in the limit A; — )• oo we have 



The main conclusion is that for A; = 4, 5, . . . , oo the fast-diffusion mean-field theory 
should be asymptotically correct for large times. Successful depositions are so far 
and few between, that the system has enough time to equilibrate diffusively between 
depositions. The point brought up in Section 3, regarding the use of the effective 
"renormalized" rate to actually fit the data by the mean-field functional forms, must 
be kept in mind, however. The case A; = 3 is borderline. Our estimates "from within 




(5.6) 



Note that we used relation (2.2) to replace p-dependence by ^-dependence. 




(A; = oo) . 



(5.7) 
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Fig. 2. The functions defined by the right-hand sides of (5.6), for k = 2,3,4,6, and 
(5.7), for k = oo. 



mean-field" cannot be used to decide if fluctuations will be marginally relevant in 
this case, i.e., if additional logarithmic factors will make the appropriately modified 
version of (5.6) break down at large times. Phenomenologically, however, the mean- 
field expressions with "effective" rate constant may provide a good quantitative fit of 
any data set taken over a large but finite number of decades in T. 
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For k = 2 the mean-field theory breaks down similarly to chemical reactions [9]. 
Other approaches which account for local fluctuations must be adopted; see [4] for 
such developments and numerical tests. 

6. DISCUSSION 

In this last section we begin by addressing the issue of how various ingredients of 
the theory are modified for d > 1. Firstly, fluctuations are always smaller in higher 
dimensions and therefore mean-field theories would provide a useful tool. However, 
as mentioned in Section 1, there are several points of difficulty. Locally, in higher 
dimensions there can be immobile "gridlocked" gaps (for nonspherical object shapes 
or on lattices), and it is not even obvious if the large time coverage will always reach 
the maximal, fully crystalline packing value. In fact, the final state may be amorphous, 
or polycrystalline with a network of defects. 

Indeed, polycrystalline and amorphous structures were found in numerical ex- 
periments on the formation of packed structures in two and three dimensions [14] by 
methods other than random sequential filling with diffusional relaxation. (We are not 
aware of any numerical investigations of the latter process in (Z > 1.) In fact, the 
tendency was [14] for polycrystallinity in 2d and amorphous structures in 3d. 

Globally, neither the dynamics of hard-core diffusing extended objects nor the 
collective aspects of the equilibrium state are well understood, even without added 
deposition [8]. Glassy slowing down and phase transitions are among the collective 
effects possible. Thus, even within the fast diffusion approximation, the problem is 
much more complicated than in one dimension. 

Another aspect of the difficulty in extending the theory to (Z > 1 is that numerical 
data are quite difficult to collect, and the currently available experimental results [6] 
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are few and very recent. However, the complexity of the problem also suggests that 
new interesting effects may be discovered. As a first step, heavy numerical effort is 
called for. 

In the present work, we accomplished the formulation of the mean-field approach 
directly for the deposition problem, without having to map to the chemical reactions. 
The main improvement was the well defined large-A; behavior of the new mean-field 
expressions, and a larger regime of validity for each value of A; > 3. The criteria of 
applicability of the mean-field theory were derived in the deposition-model context. A 
new finding was that mean-field theory provides the correct asymptotic description of 
continuum deposition, with diffusional relaxation, on a line. 
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